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Abstract. - We show that a 3-qubit system as studied for quantum information purposes 
can alternatively be used as a thermodynamic machine when externally driven in finite time 

"t^ ' and interfaced between two split baths. The spins are arranged in a chain where the working 

spin in the middle exercises Carnot cycles the area of which defines the exchanged work. The 
cycle orientation (sign of the exchanged work) flips as the difference of bath temperatures goes 

i-rt ' through a critical value. 

o 

Back in 1824 Sadi Carnot published his "Reflections on the Motive Power of Fire" , propos- 
CN ■ ing ideas for a combustion engine [1]. His reasoning was based on his now famous ideal gas 
cycle analysis. With the advent of nano physics the intriguing question arises: Down to which 
scales can one reasonable define thermodynamic variables and thermodynamic machines? 
^^ I First attempts on the quantum level have been made in [2,3]. Meanwhile such miniature 

^^ ■ machines have attracted quite some attention [4-9]. Artificial autonomous nanomotors pow- 
^^ I ered by visible light have recently been demonstrated experimentally [10]. Nanomechanical 
devices have been of growing interest in the last couple of years [11]; the importance of cooling 
processes is being recognized with respect to quantum information processing [12]. Further- 
jrt ' more, it has been realized that strange analogies appear to exist between the basic concepts 
of quantum information theory like specific entanglement measures on the one hand side and 
thermodynamic concepts on the other side [13], the origin of which is still fairly unclear. 
' O I Two level systems (TLS) like spins or qubits are the essential ingredient for quantum 

computation. Much effort has been directed towards control of small clusters and chains 
of qubits in quantum optical systems [14], nuclear magnetic resonance [15] and solid state 
systems [16]. A serious problem in any such realisation is the interaction of the respective 
quantum network with its environment. 

Here we argue that structurally identical quantum networks can be used as quantum 



JH I information processors as well as quantum thermodynamic engines! This may not only shed 
new light on that intricate relationship, but may also give new directions for the application 
of such few-TLS structures, for which a number of different implementations are becoming 
available now. 

(*) E-mail: Markus.Henrichaitpl.uni-stuttgart .de 
© EDP Sciences 
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Fig. 1 - Schematic representation of the system under investigation: An inhomogeneous 3-spin-chain 
is interfaced between two baths while a time-dependent field deforms the spectrum of spin 2. 



The model under consideration will be a Heisenberg coupled spin-chain consisting of three 
spins which is a typical theoretical description for TLS. Its Hamiltonian reads 
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(1) 



where fi specifies the spin, Sf^ its local energy splitting, A the nearest neighbour coupling 
strength and CTi(/z) are the Pauli-matrices. The detuning between spin 1 and spin 3 is S13 = 
{Si — 5z)/'2 > 0. This model was chosen because of its analogy to the classical carnot engine. 
As will be explained below the spin in the middle can be considered as the "working gas" 
whereas the other two spins act as contacts and filters. These filters introduce a resonance 
effect which is essential for the model to work as a kind of carnot engine. 

Thermodynamic properties can be imparted on this system by an appropriate embedding 
into a larger quantum environment [17-19], without the need of any thermal bath. In the 
present context, though, it is much simpler to settle for the customary approach: As depicted 
in Fig. [1] the chain is taken to be in thermal contact with two baths a = c,h with different 
temperatures Tc < Th, such that spin 1 is weakly coupled to the hot bath only and spin 3 to 
the cold bath. A convenient method to simulate such an open system is then based on the 
quantum master equation. 

To derive the master equation for our model one starts from the Liouville-von-Neumann 
equation (we set h = 1 and also the Boltzmann constant /cg = 1). A typical bosonic bath 
model consists of uncoupled harmonic oscillators. With the help of projection operator tech- 
nique up to second order and the Born-Markov approximation [20] the master-equation for 
the whole spin system then takes the form [21] 



dt 



p{t) = -i[H, p{t)] + Mm) + -Dcim) 



(2) 



with H according to ([1} and p denoting the density matrix of the spin system. The dissipator 
I?Q for bath a = c,h is 



+(p(i)f ii„) ~ (i„/5(i)f 1) , 

where Aa are the local coupling operators 

A = CT,(i)0i(2)®i(3) 



(3) 



(4) 
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i, = 1(1)® 1(2)® (7,(3). (5) 

The matrix elements of Fq, in the eigenbasis have the form [18, 21] 

(^1 r« 1^-) - - (^5?^ + ^('^-)^S^) ^S- (6) 

Here k is a couphng parameter between system and either bath a, A" are the matrix elements 
of Aa, 0{uj) is the step function, /3q, = 1/Ta the bath temperature and ujij = Ei — Ej are 
the frequencies of the total spin system {H \i) = Ei \i)). The stationary state of this master 
equation with only one bath and for fixed i5^ is easily shown to be canonical of the form 
/3s = e~^^ Tr {e^''^}. With both baths in place the spin chain might be viewed as a molecular 
bridge generating a stationary leakage current Jl = Jh = —Jc- 

Here the heat current J^ between the 3-spin-systcm and the bath a can be defined by the 
energy dissipated via bath a [22]: 

J^ = Tr {HV^ipit))} (7) 

where Tr means the trace operation. In the following a current out of bath a into the machine 
will be defined as positive. 

This situation changes drastically, if in addition, the center spin 2 is driven externally: We 
assume the form S2{t) = sin{ujt)b + bg with oj ^ (5^, which enables the bath to sufficiently 
damp spin 2. If spin 2 was driven too fast it would no longer be possible to exchange energy 
between spin 2 via spin 1 (3) and its respective bath. The parameters are chosen such that 
(52 (t) oscillates between the two resonance conditions, i. e. bo = {Si + 5^)/2 and b is the 
detuning (^13. It should be possible to implement such devices with any TLS which allows 
energy exchange based on a coupling small compared with the local energy splittings. In case 
of an actual spin (52(t) could be implemented by a varying magnetic field B. If the TLS was 
realized by a semiconductor quantum dot, the spectrum might be distorted by an external 
electric field E. 

For any diagonal state of a TLS fi in its local energy basis we can always define the tem- 
perature T^ = —Sfi/ [in {pji/pfi)] and the von Neumann entropy S'^ = —{pf^ lnp° + pj^ InpjJ, 
where p° are occupation probabilities for the ground state and p^ for the excited state. Under 
the conditions as stated there can be no local non-equilibrium, a restrictive but convenient 
property of TLS. As will become transparent below, T2 and S2 indeed behave like thermody- 
namic variables. 

Now we know from the Gibbs relation for the energy change of subsystem 2 that [23] 

AE2^AQ2 + AW2^-J2Ap^S^2 + ^T.P^^^'^- (8) 

The heat AQ2 is thus associated with the change of occupation probabilities P2, the work 
AW2 with the deformation of the spectrum by external driving (modulation of 52{t)). Over 
one cycle AE2 = so that 

W2 = -Q2 = - i T2dS2. (9) 



On the other hand the energy fiowing from or into bath a is heat by definition. Over one 
period t — I-k juj this heat is 

/•T 

J„dt (10) 
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Fig. 2 - ^Tb-Diagram for the quantum heat pump. The arrows indicate the direction of circulation. 



allowing us to identify Q2 = Qc + Qh- The work per cycle, W2, put into the machine will be 
counted as positive. This sign convention also applies to the transferred heat per cycle, Qa- 

All numerical results are obtained by integrating the master equation with a 4th order 
Runge-Kutta algorithm. Because of (52 (i), H is also time-depended. Thus for each time step 
one has to calculate new dissipators 2?q. The parameters taken here are A = 0.01, k = 0.001, 
^1 = 2.25, ^3 = 1.75, w = 1/128, Tc = 0.4 and Th is varied, unless stated otherwise. Both 
coupling parameters A and n are chosen to stay in the weak coupling limit. Figure [5] shows 
the dynamics of subsystem 2 in its entropy-temperature space. 

In analogy with the classical Carnot machine one may approximately distinguish four 
different steps in the S'2T2-diagram (Fig. [2]) in combination with the respective currents in 
Fig. El 

1. Quasi-iso-thermal step: Spin 2 (the "working gas") is brought in resonance with spin 
3 {S2{t) « 63) and so couples with bath c at temperature Tc. Because of this energy 
resonance the current Jc via spin 3 will increase, whereas the current Jh via spin 2 will 
be negligible. The occupation probabilities of spin 2 and 3 will approach each other and 
so do the respective local temperatures. 

2. Quasi-adiabatic step: Spin 2 is driven out of resonance with spin 3 {Si > S2{t) > (S3), 
now Jc decrease while Jh nearly stays unchanged. The occupation probability of spin 2 
does not change significantly and there is almost no change in the entropy 5*2. 
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Fig. 3 - The heat current Ja(i) over one cycle with duration r = 27r/a; — 804.25. The peaks are a 
result of the resonance-effect. 
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Fig. 4 - The Carnot-efRciency 17^3^^ and the efficiency rj^j^ of the quantum heat pump (AT < ATmax) 
and rjQ^,. and rj^^ of the heat engine ( AT > ATmax) as function of the temperature difference AT. 
'7ond — 1 ^ \/Tc/Th is the efficiency of the classical endoreversible heat engine [24] . 



3. Quasi-iso-thermal step: Spin 2 comes in resonance with spin 1 (<52(t) ~ Si) and by that 
in contact with bath h at temperature T^. Jh increases whereas Jc is very small. The 
local temperatures of spin 1 and 2 nearly equal each other. 

4. Quasi-adiabatic step: Spin 2 is driven back to the starting point as in step 2. 

Also under driving some leakage current Jl will flow from the hotter bath h to the colder 
bath c leading to a transferred heat Ql — § Jh^t per cycle. This heat is always contained 
in Qc, Qh with the proper sign. W2 = —{Qh + Qc) is not influenced by Ql- 

In contrast to the ideal Carnot heat pump which is defined in the quasi-static limit and 
with the bath coupling set exactly to zero during the adiabatic steps, the quantum heat 
pump must be driven through its cycle in a finite time (cf. the model of the endoreversible 
engine [24] ) . If one would wait after each infinitesimal time step till the system had reached its 
momentary steady transport configuration (i. e. w <C k), the quantum heat pump would no 
longer work. The area of the S'2T'2-cycle would be zero. In step 1 (or 4) the local temperature 
of spin 2 is larger (smaller) than the respective bath temperature Th (Tc) which is needed to 
drive a current from spin 2 to bath h (from bath c to spin 2). There is an entropy production 
due the dissipation Va- 

The efficiency of a heat pump is defined by the ratio of the heat Qh pumped to the hot 
reservoir and the work applied 77^ — Qh/W2,'which. reduces for the Carnot heat pump to 
^Car ~ T^h/iTh — Tc). For the heat engine the efficiency is 77^ = W2/Q/1 which in the Carnot 
case leads to 775^^,. = 1 — Tc/Th- 

For the quantum heat pump we have calculated the work 11^2 for one cycle out of the 
S'2T2-diagram by (H]), and the heat Qa by (HO]). We find that indeed VF2 + Qc + Qh = in all 
cases and by that confirm the use of T2 and S2 as effective thermodynamic variables. 

Figure m (left part) shows the numerical results for the efficiency 77^^ as function of AT = 
Th — Tc together with the corresponding Carnot efficiency 77^^^. . As can be seen, the efficiency 
of the quantum heat pump is always smaller than that of the Carnot pump, as should have 
been expected. Note that for AT ~ the efficiency does not diverge as the Carnot efficiency 
does. Enhanced cooling extends to AT < 0, i. e. Jc can be made larger than the reversed 
leakage current. 

The work W2 to be applied goes to zero at AT = ATmax (cf- Fig. [S]). A rough estimate 
gives ATmax = Tc {61/63 — 1) Ri 0.71 for our model parameters. The actual numerical value 
is 0.60. At this point the mechanical driving becomes lossless, only the leakage remains. The 
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Fig. 5 - Heat Qc and Qh as well as the work W2 performed over one cycle as function of the temper- 
ature difference AT. The inset shows these functions around the point AT = ATmax in more detail. 
Ql is the leakage heat per cycle. 



efficiency 77^^ goes to zero at some AT < ATmax, when the heat pumped up is just cancelled 
by the heat leaked down, i. e. when Qh = 0. 

To make this conclusion clearer, also the heat Qc and Qh per cycle are included in Fig. 
[S] as function of the temperature difference. The inset shows the point AT — ATmax for 
which W2 = and Qh = Ql = —Qc- Qc changes its sign before Qh and Qh before W2 as 
AT approaches ATmax from below. For AT > ATmax the machine starts working as a heat 
engine {W2 < 0) with efficiencies rj^^ ('7car) shown in Fig. |4] at right. ?7qm is not far away 
from the classical efficiency of the so-called endoreversible heat engine, i. e. for a machine 
completing its cycles at finite time [24] . Eventually this machine operates irreversibly, as the 
external heat flow from or into the baths is modelled as a transport process. Figure [6] shows 
the corresponding S'2T2-diagram, now executed with the opposite orientation, in agreement 
with the reversed sign of W2 ■ 

Before concluding we remark that these 3-qubit-systems could also be used as nano-scale 
sensors: If ATmax was known as a device property, the temperature difference AT between 
an unknown bath Tc and a know bath Th could be determined by varying Th until W2 = 0. 
At this point the measurement was least disturbing. 

Driven spin networks have long since been basis for elementary quantum information 
processing tasks. For these the influence of decoherence (baths) is disastrous. In this letter 
we have argued that, when (two conflicting) thermal contacts are added deliberately, such 
networks should show a completely different, but still quite fascinating behavior. This might 
open up a new field of research: general quantum networks under the influence of unitary gates 
(interface to work reservoirs) and thermal contacts (interface to heat reservoirs). In a further 




Fig. 6 - S'2r2-Diagram for the quantum heat engine. The arrows indicate the direction of circulation. 
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twist also those reservoirs may be taken to be finite and described by quantum mechanics in 
full detail. 
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